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DISCLAIMER 



Although each program has been tested by its contributor, 
no warranty, express or implied, is made by the contributor 
'or any User^s Group, as to the accuracy and functioning of 
the program and related program material, nor shall the fact 
of distribution constitute any such warranty, and no responsibility 
is assumed by the contributor or any User's Group, in connection 
therewith. 



COMMON USER3 GROUP PROGRAM REVIEW AND EVALUATION 
(fill out in typewriter, ink or pencil) 

Program No. . Date _^ 

Program Name: 

1. Does the abstract adequately describe what the program is and what Yes No_ 

it does ? 

C om m ent 

2. Does the program do what the abstract says ? Yes No_ 

Comnient ^ 

3. Is the description clear, understandable, and adequate? Yes No_ 

Comment 

4. Are the Operating Instructions understandable and in sufficient detail? Yes No^ 

Comment ^ 

Are the Sense Switch options adequately described (if applicable) ? Yes No_ 

Are the mnemonic labels identified or sufficiently understandable ? Yes No_ 

Comment 



5. Does the source program compile satisfactorily (if applicable)? Yes No_ 

Comment 

6. Does the object program run satisfactorily? Yes No_ 

Comment ■ 

7. Number of test cases run . Are any restrictions as to data, 

size, range, etc. covered adequately in description? Yes No_ 

Comment . 

8. Does the Program meet the minimal standards of COMMON? Yes No_ 

Comment 

9. Were all necessary parts of the program received? Yes No 

Comment__ 

10. Please list on the back any suggestions to improve the usefulness of the program. 
These will be passed onto the author for his consideration. 

Please return to: Your Name 

I3M Corporation Company 

Program Information Department • ■ . 

40 Saw Mill River Road Address 



Hawthorne, New York 10532 ^ - ■ 

Attn: PREP FORM COORDINATOR Users Group Code 



THiS REVIEW FORM IS PART OF THE COMMON ORGANIZATION'S PROGRAM REVIEW ANI 

EVALUATION PROCEDURE. NONMEMBERS ARE CORDIALLY INVITED TO PARTICIPATE 
IN THIS EVA LUATION. 
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Factor Analysis by the Centrold 

Technique 
William F., Mllnei- 
North Texas State University 
Kay 10, 1962 
6,0 statistical Program 



Modifications or revisions to this program, as they occur, 
will be announced in the appropriate Catalog of Programs 
for IBM Data Processing Systems. When such an announce- 
ment occurs, users should order a complete new program 
from the Program Information Department. 
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il. Detailed Discrlption cf rrccran: 

III. Input / Output 

IV. Senple Problem 

V, ^pdratlng Instructions 

VI, Program Listing 



FvoRTam Hanxial 

I. Title: Factor Analysis by Centrold-NTeoeS 

II. Description of Program 
A. aes\ilts 

The program output ia the factor loadings associated 
with the raw data. A first approximation is typed 
first and then as a second pass a refined set of 
factor loadinfrs, based on the first approxinatlons, 
are output. 
3, Kethod of Crsputalion 

Itiurston's method cf ccntrcM 

III. Input/Cutput 
Input Por?iat 

Card 1: il blank iC2 

where fH = the number of observed variables 

H2 = the number of obs6z>vations 
Card 2: R^j ,1=1, 3=1 
Card 3i R12 



Output Format 

The factor loadings matrix in row order 

The as described by Thurston in row order. 



n 

.00000000 fn 

,0000005 fix 

,59229011 r,3 

.266if7283 ri •/ 

1.0000005 ^, Stmfl^ (MiciihAJ 

.00000000 f:,^ Cofff^cK^rMfny 

,59228991 iH ^UJ drdif 

,26647288 6»/ it) it\\ IttXJ ^ ^^Q*^^ 

,59229011 1^, 

,59228991 

.00000000 

,501»32«*03 Kji^ 

.26647283 K^, 

.26647288 

,50432403 

.00000000 
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5.905'»'»85E-02 
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-.W»751867 

.3923013'* 
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-5;<^5281»8E-02 
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1;0%3<»253 
1,0508116 
-.66570587 
.*7 137924 
STOP 



? 

• Operating Instructions 

1) Set tabs and maragln as for standard POBtTRAN 

2) Ready card reader with program deck followed Ay 
data cards 

3) Depress reset and load, when the "Load Data" 
awssage is typed depress start on the 1620 

{{.) Tha program will read the first data card 
and atop ready to accept two p4ice8 of data 
fi*aiii the typewriter. These are 

KAP the number of factors to extract up to 12 
Data the amount of comunallty to extract j this 
•la rarely known and a value of 1000, is 
adequate to estract all. 
5) Tto pro/s^ram will then read the balance of the data 
cai>ds and will output the first approximation 
for the factor loadings. After typing these the 
program will attempt to read again. If a refined 
set of loadings and the are desired place the 
entire data deck In the reader and hit start. 



C FACTOR ANALYSIS, CENTROID 

C INPUT FOR THIS PGM IS OUTPUT OF THE SIMPLE CORREUTION 

C ROUTINE BY THE SAME AUTHOR. ALL DATA CARDS WILL BE 

C AND THE PGM WILL COMPUTE AND OUTPUT THE FIRST APPROX. 

C THE FACTOR- LOADINGS, THE PGM WILL THEN TRY TO READ THE 

C DATA AGAIN. AT THIS TIME THE ORIGINAL DATA DECK SHOULD 0^, . . 

C ENTERED AGAIN AND THE OUTPUT WILL BE A REFINED SET OF /^c^^i^MOj- 

C GENE MILNER 

C NORTH TEXAS STATE UNIVERSITY 

DIMENSION R{I5,15),C{15,12),SUMI(15),H2(15) 

1 REA0,N1,N2 
ACCEPT, KAP .DATA 

2 00 3 1-1, n! 

3 H2(l)-0. 
MAP-1 
MIG*1 
H-0. 
MUG-I - 
SUMH-0. 

DO k i-1,N1 00 k l>1,N1 

DO k J«1,N1 
k REA0,R(l,J) 
500 K-1 

00 510 l-I.NI 

DO 510 Jiil,12 
510 C(I,J)-1.0 
525 L-I 

530 00 5 1-1. N1 

5 SUMI(n-0. 
SUMC-0.- 
SUMC2-0. 

DO 6 l-1,N1 
DO 6 J-1,N1 

6 SUN1(I}»SUMI(I)+R(I,J) 
T-O. 

00 J l-l,N1 

7 T-t+SUMI(l) 
N4-N1+1 

IF(L-1)725,590,725 
590 00 8 1-1. n! 

IF (SUM! (I))610,8,8 

8 CONTINUE 
GO TO 670 

610 6A-990«». 
DO 9 1-1, Nt» 

IF(SUM|{I) - GA) 620,9,9 
620 GA-SUMI(I) 
JIG-I 

9 CONTINUE 
l-JIG 

00 10 J-1,N1 

10 R(I,J)— R(I,J) 
J-JIG 

00 11 l-1,N1 

11 R(I,J)— RU.J) 
C(J|- ' ■ 



ltG,K)— C(JIG,K) 
GO TO 530 
670 60 TO (67 1.701), MUG 



o 



671 DO 700 1-1, N1 
RAG-0. 

DO 690 J-1,N1 

IF(ABS(R(I,J))-RAG)690,690,685 
685 RAG-ABS{R(I,J)) 
690 CONTINUE 

700 R(I,I)-RAG 

701 MUG-1 

GO T0(705,72O),MAP 
705 MAP-2 

00 715 l«1,N1 
715 H-H+R{l,l) 
720 L-2 

60 TO 530 
725 IF(L-2)730. 770,730 
730 DO 7kO 1-1, Nl 
735 IF(SUMI{I))750,7'»0,750 
7'»0 CONTINUE 

60 TO 755 
750 CONTINUE 

60 TO 7'»0 
755 00 765 1*1, Nl 
765 R(l,l)=0. 

GO TO 525 
770 RTT-SQR(T) 

DO 805 1-1, Nl 

TEMP-SUMI{I)/RTT 

IF(C(I.K))785,790,790 
785 IF(TEMP)795,800,800 
790 IF(TEMP)800,795,795 
795 C{I,K)-TEMP 

GO TO 805 
800 C(I,K)— TEMP 
805 CONTINUE 

DO 815 l-1,N1 

SW1C-SUMC+ABS(C(I,K)) 
815 SUMC2-SUMC2+C(I,K)*C(I,K) 

IF{SUMC-RTT)820,825,820 
820 TYPE,SUMC,RTT 
825 PHK-SUMC2/H 

Sl*1H-SUMH+PHK 

TYPE,PHK 

DO 12 1-1, Nl 

TYPE,C(I,K),I,K 

12 CONTINUE 

DO 835 l-1,N1 
00 835 J«1,N1 
835 R(i,J)-R(l,J)-ABS(C(l,K)*C(J,K)) 

I F ( S WH-OATA) 81»0 ,8»» . 8 5 5 
8I»0 IF(KAP-K)855,855,850 
850 I&-K41 

60 TO 530 

855 GO TO(856,900),MIG 

856 00 865 1-1. Nl 
860 DO 865 J«1,K 

865 H2(I)-H2(I)+C(I,J)*C{I,J) 
READ,N1,N2 
DO 15 1-1, Nl 
00 13 J-1,N1 
REAO,R(i,J) 

13 CONTINUE 



DO 885 1-1, Nl 
885 R(I.I)-H2(I) 

HUG-2 

MAP*1 

H-0. - 

SUMH-0. 

Ml 6-2 

60 TO 500 
900 CONTINUE 

DO 16 I-1,N1 
^ TYPEjH2(lJ 
16 CONTINUE 

STOP 

END 



